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Abstract. We introduce a new method to achieve long lived Bloch oscillations and 
dynamical localization of matter wave gap solitons in optical lattices. The method is 
based on time dependent modulations of the nonlinearity which can be experimentally 
implemented by means of the Feshbach resonance technique. In particular, we show 
that the width of the wave packet is preserved if time modulations of the nonlinearity 
are taken proportional to the curvature of the linear band spectrum which for most 
typical experimental settings are well approximated by harmonic time modulations of 
proper frequencies. 
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Existence of non spreading localized waves (or solitons) is a property of nonlinear 
systems which steams from the balance of nonlinearity and dispersion. In the case of 
Bose-Einstein condensates (BECs) tightly bounded in two spatial dimensions, bright 
matter solitons are known to exist when attractive two-body interactions compensate 
the vacuum dispersion. These waves were experimentally reported in Ref. pQ. An 
external periodic potential generated by two counter propagating laser beams, i.e. an 
optical lattice (OL), applied to a BEC gives rise to a band-gap structure in the spectrum, 
introducing in this way artificial dispersion. In the effective mass approximation [2] one 
finds that signs of the atomic effective masses are different at opposite edges of each finite 
gap. Under this condition the properties of matter waves change significantly: a BEC 
with repulsive inter-atomic interactions admits the existence of stationary gap solitons 
with chemical potentials inside spectrum band gaps. The experimental observation 
of these matter wave gap solitons was reported in [3]. Due to broken translational 
symmetry these solutions are intrinsically localized, especially for high strength OLs, so 
that any attempt to move them with the aid of an external force (such as, for example, 
the acceleration of the OL), usually leads to their destruction. 

In order to control gap solitons and allow them to move in real and in reciprocal 
spaces, one has to properly design the nonlinearity in the system . Soliton existence 
requires indeed that the effective nonlinearity and the effective mass have opposite signs 
for all points in the Brillouin zone [2]. In recent papers [U [5] it has been shown that 
this condition can be satisfied with the help of a nonlinear optical lattice, i.e. with a 
spatial modulation of the scattering length, and that it leads to the existence of Bloch 
oscillations and dynamical localization of matter waves which are long lived in presence 
of nonlinearity. A nonlinear optical lattice imply a control of the scattering length in 
space, a task which is more difficult to implement than the modulation of the scattering 
length in time (this last can be realized with external time dependent magnetic fields 
via the usual Feshbach resonance technique [6]). It is then natural to ask if a time 
management of the nonlinearity can lead to the above long lived Bloch and dynamical 
oscillatory behaviors in the nonlinear regime. 

The aim of the present Letter is just devoted to this. In particular we show 
for the first time, that long lived Bloch oscillations and dynamical localization of a 
matter wave gap soliton in an accelerated linear OL become possible in presence of a 
properly designed time management of the nonlinearity. We remark that simultaneous 
effects of linear OLs and nonlinear time-dependent nonlinearities has been considered 
in literature in various contexts. In the case of BECs this technique was used mainly 
for the stabilization of higher dimensional solitons [7], investigation of stable Feshbach 
resonance-managed discrete matter-wave solitons [8] and generation of bright and dark 
solitons [9]. No applications of this technique to nonlinear Bloch oscillations and 
dynamical localization phenomena has been suggested so far. We remark that the 
implementation of properly designed time dependent nonlinearities for non destructive 
matter wave transport in OLs can be of interest non only for BECs but also for nonlinear 
optics. 
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Let us start by considering a BEC described in the mean-field approximation by the 
one-dimensional Gross-Pitaevski equation (GPE) with a time-varying s-wave scattering 
length a s (T) and cos-like OL accelerated in space according to the law S(T): 



lh df = -2^dX- 2 - Vc ° S 



^(X-S(T)) 



^ + ^^l^l 2 ^. (1) 



a ± m 

Here m and \P denote the boson mass and the macroscopic wavefunction, V is the OL 
amplitude, T and X are the time and the longitudinal coordinate, respectively, d is 
the OL period, a± is the transverse trap size. Eq.([T]) is obtained from the full three- 
dimensional GPE in the case, when OL period d is much bigger than transversal trap 
size aiCd (for details of transition see, e.g., [9]). Introducing the dimensionless moving 
spatial coordinate x = ir(X — S(T))/d, time t = TE R /h (E R = h 2 ii 2 / {2md 2 ) is the recoil 
energy) and wavefunction 



ip{x,t) = — W(a, Tjexp 



.mdS ,,_ - m f T ( dS\ 2 , 

,__(jr-s ( r)) + ,_jf dx 



we obtain the normalized GPE 

= ~ 1 TT-vcos{2x)iP + 'y{t)x^ + g{t)\i;\ 2 ^ (2) 



with dimensionless OL amplitude w = V/Er, time-dependent nonlinearity g(t) = 
a s (T)/\a s (0)\ and external force j(t) = (md/irE R )d 2 S/dT 2 , proportional to the 
acceleration of the OL. In the absence of the nonlinearity g(t) = and external force 
7(t) = 0, periodicity of OL gives rise to band-gap structure of both spectrum e n (q) and 
Bloch functions <p nq (x) of the linear problem e n (q)ip nq = —d 2 ip nq /dx 2 — v cos(2x)ip nq , 
where n is the band number ( below we consider only the first band n = 1 and therefore 
the index n is omitted) and q is Bloch wavenumber. The spectrum is periodic in the 
reciprocal space with the period 2 (in the chosen units): e(q) = e(q + 2). 

Under the influence of the external force, soliton motion in the real and reciprocal 
spaces obeys semiclassical laws for its center in real and reciprocal space X and Q, 
correspondingly [U [5] : 



X = — 
dq 



, Q= -7(0, (3) 

q=Q 

where the overdot stands for the time derivative. When the external force causes soliton 
motion, as it follows from ([3]), the position of soliton center Q in reciprocal space 
will be changed. At the same time the existence of a gap soliton solution of Eq.Q 
for a given Q{t) = Q — J*j(r)dr (Qo is the initial coordinate of the soliton in the 
reciprocal space) is determined by the condition M(Q(t))g(t) < for all times t (here 
M(q) = [d 2 e(q)/dq 2 ] 1 is the soliton effective mass). Thus, to prevent soliton destruction 
one should vary the nonlinearity g(t) in a manner to keep opposite sign with respect the 
effective mass, (this assures the above existence condition is satisfied for all t and all 
Q{t)). Moreover, to facilitate the soliton motion, one should minimize redistribution of 
particles, keeping the soliton width to be constant (this can be achieved by keeping the 
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Figure 1. (Color online) (a,d) Dependencies g(t), providing condition ^ for OL 
amplitude v = 3, external force T = —0.001, and initial Bloch wavenumbcrs Q = 1 
(panel a), Qq = (panel d). (b,e) Long-living Bloch oscillations in BEC with time- 
varying nonlincarity, depicted in panels (a) and (d), respectively. (c,f) Decaying Bloch 
oscillations in BEC with the same parameters, as that of panels b,e, respectively, 
except constant nonlincarity, which is g(t) = 1 in (c) and g(t) = —1 in (f). The 
initial condition is the stationary gap soliton with E = —0.7324, Qq = 1 (b,c) and 
E = —0.938, Qq — (c,f). All parameters are chosen the same as in @]. 



product of effective mass and effective nonlinearity g(t) J_ \ipQ^(x)\ A dx to be constant 
jlj). For this it is sufficient to take g(t) of the form 

= \M(Q )\j: n \y Qo (x)\ A dx 

9[) Mm)) /:>Q (i )(*)iw u 

In the following we apply the above nonlinear management scheme to induce long lived 
nonlinear Bloch oscillations and dynamical localization of a gap soliton in an accelerated 
OL. 

Bloch oscillations of matter waves occur when BEC is subject to a constant (dc) 
external force, i.e. for j(t) = T =const. For the particular choice of the lattice amplitude 
and the linear force, the above algorithm of designing the temporal dependence (jl]) of the 
nonlinearity results in curves depicted in Fig. [T^,d for Q — 1 and Qq = 0, respectively. 
The gap-soliton dynamics under this g(t) management is depicted in Fig. [Tb,e. Notice, 
that due to proper choice of g(t) Bloch oscillations become very regular with the period 
and the spatial amplitude which coincide with the values T s = 2/\T\ = 2 • 10 3 and 
X s = [e(l) — e(0)]/(2|r|) ~ 32.57T predicted from Eqs.([3]). This behavior is in deep 
contrast with the one depicted in Fig. [Tb,f for the case of a constant nonlinearity, where 
the Bloch oscillations are accompanied by the decay of the wave packet. 

Similar results are obtained for the case of the dynamical localization. In this case 
the motion of a soliton is turned on under the action of a time-periodic external force 
7(i) = r cos (Qt). When the frequency Q is much smaller than the soliton energy E 
(approximated by the chemical potentials of the initial stationary state for the initial 



Figure 2. (Color online) Dependencies g(t), providing condition ^ (panels a,c) and 
correspondent particle densities |^(x,f)| 2 (panels b,d, respectively) for OL amplitude 
v = 3, initial Bloch wavenumbcr Qo = 1 and external force with parameters 
r = -0.001, t = 500 and frequencies w 0.001306 (panel a,b) and SI w 0.000569 
(panel c,d). All parameters are chosen in accordance with [5]. 



Qo = 0,1) and at the same time fulfills the condition Q = Tir/z n with z n being the n-th 
zero of the Bessel function Jo(z n ) = [10], one expects dynamical localization of the 
soliton, i.e. a regular oscillation in space induced by the ac force, provided its existence 
is guaranteed for all the energies belonging to the chose band. In [5] this was achieved 
with help of spatially non uniform nonlinearity. Now we design the time-dependent 
nonlinearity in order to obtain the same effect. 

Application of external force, which starts to oscillate at time t — to 



and time-dependent nonlinearity g(t) (as depicted in FigfSk^c) results in the dynamical 
localization of the soliton, when its mean drift velocity becomes to be zero. This results 
are obtained both for the oscillation frequency, which corresponds to the first zero of 
Bessel function (Figj2b) and for oscillation frequency, determined from the second zero 
of Bessel function (FigfSH). 

The numerical simulations, reported above, are performed for exact numerical 
design of g(t). However for typical values of v > 1 they are remarkably well reproduced 
using the dependence g(t) which steams form the tight-binding approximation. Since in 
typical experimental settings v is of the order of several recoil energies, expanding e(q) in 
the Fourier series, we can consider only the two leading terms: e(q) = u — oj\ cos(7rg). 
This approximation works well already at v > 1, as one can see form the Table I of 
Ref. pj]. For the particular case v — 3 (for which our above-mentioned numerical 
simulations are performed), c^o ~ —0.84, oj\ ~ 0.102. On the other hand, in this 
approximation the integral in Eq. PJ) becomes practically a constant, independent on 
q and, hence, on £0 and the temporal law (j4j) for nonlinearity g(t) simplifies to 



f The deviation from a constant of the integral in ((4]), due to its dpendence on Q, can be estimated 
in the Wannier basis as the overlapping integral of Wannicr functions of two adjacent sites, this giving 
for v = 3 an error of about 0.2 %. 




(5) 




(6) 



6 

i.e. it is uniquely fixed by the dynamics of the wavepacket in the reciprocal space. 

In this case of Bloch oscillations the semiclassical equations of motion give X(t) = 
X — jr cos(nQo) [1 — cos(7iTt)] and Q(t) = Qo — Tt with X the initial coordinate of 
the soliton. The temporal dependence of the nonlinearity function then follows from 
Eq. (EI) as g(t) m — cos(ttQo) cos(7rrt). At the same time for the dynamical localization 
case time-dependent nonlinearity can be approximated from Eq.Q as 

g(t) „ - co S (,rQ ) j ^ ^ + T _ ^ ^ _ ^ ( £ ^ . (7) 

To conclude, we have introduced a new method based on time dependent 
modulations of the nonlinearity to achieve long lived Bloch oscillation and dynamical 
localization of gap solitons in OLs. The desired modulation is obtained directly from 
the curvature of the underlying band structure and can be experimentally implemented 
by means of time dependent external magnetic fields via the Feshbach resonance 
technique. For strengths of the OL greather than one recoil energy (e.g. for most 
typical experimental settings), the shape of the time modulation can be fixed in very 
simple manner using fields of proper frequencies. 
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